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Abstract — New upper and lower bounds are presented on the 
capacity of the free-space optical intensity channel. This channel 
is characterized by inputs that are nonnegative (representing the 
transmitted optical intensity) and by outputs that are corrupted 
by additive white Gaussian noise (because in free space the distur- 
bances arise from many independent sources). Due to battery and 
safety reasons the inputs are simultaneously constrained in both 
their average and peak power. For a fixed ratio of the average 
power to the peak power the difference between the upper and 
the lower bounds tends to zero as the average power tends to 
infinity, and the ratio of the upper and lower bounds tends to 
one as the average power tends to zero. 

The case where only an average-power constraint is imposed 
on the input is treated separately. In this case, the difference of 
the upper and lower bound tends to as the average power tends 
to infinity, and their ratio tends to a constant as the power tends 
to zero. 

I. Introduction 

We consider a channel model for short optical communi- 
cations in free space such as the communication between a 
remote control and a TV. We assume a channel model based 
on intensity modulation where the signal is modulated onto 
the optical intensity of the emitted light. Thus, the channel 
input is proportional to the light intensity and is therefore 
nonnegative. We further assume that the receiver directly 
measures the incident optical intensity of the incoming signal, 
i.e., it produces an electrical current at the output which is 
proportional to the detected intensity. Since in ambient light 
conditions the received signal is disturbed by a high number of 
independent sources, we model the noise as Gaussian. More- 
over, we assume that the line-of-sight component is dominant 
and ignore any effects due to multiple-path propagation like 
fading or inter-symbol interference. 

Optical communication is restricted not only by battery 
power, but also for safety reasons by the maximum al- 
lowed peak power. We therefore assume simultaneously two 
constraints: an average-power constraint £ and a maximum 
allowed peak power A. The situation where only a peak-power 
constraint is imposed corresponds to £ = A. The case of only 
an average-power constraint is treated separately. 

In this work we study the channel capacity of such an 
optical communication channel and present new upper and 
lower bounds. The maximum gap between upper and lower 
bound never exceeds 1 nat when the ratio of the average- 
power constraint to the peak-power constraint is larger than 



0. 03 or when only the average power is constrained but not 
the peak power. Asymptotically when the available average 
and peak power tend to infinity with their ratio held fixed, the 
upper and lower bounds coincide, i.e., their difference tends 
to 0. We also present the asymptotic behavior of the channel 
capacity in the limit when the power tends to 0. 

The channel model has been studied before, e.g., in [1] and 
is described in detail in the following. The received signal is 
corrupted by additive noise due to strong ambient light that 
causes high-intensity shot noise in the electrical output signal. 
In a first approximation this shot noise can be assumed to be 
independent of the signal itself, and since the noise is caused 
by many independent sources it is reasonable to model it as an 
independent and identically distributed (IID) Gaussian process. 
Also, without loss of generality we assume the noise to be 
zero-mean, since the receiver can always subtract or add any 
constant signal. 

Hence, the channel output Yk at time k, modeling a sample 
of the electrical output signal, is given by 

Y k = x k + Z k , (1) 

where x k G denotes the time-fc channel input and 

represents a sample of the electrical input current that is 
proportional to the optical intensity and therefore nonnegative, 
and where the random process {Z k } modeling the additive 
noise is given by 

{Z k }~HDN m (0,a 2 ). (2) 

It is important to note that, unlike the input x k , the output Y k 
may be negative since the noise introduced at the receiver can 
be negative. 

Since the optical intensity is proportional to the optical 
power, in such a system the instantaneous optical power is 
proportional to the electrical input current [2]. This is in 
contrast to radio communication where usually the power is 
proportional to the square of the input current. Therefore, in 
addition to the implicit nonnegativity constraint on the input, 

X k > 0, (3) 

we assume constraints both on the peak and the average power, 

1. e., 

Pr[X k > A] = 0, (4) 



E[X fc ] <£. 



(5) 



We shall denote the ratio between the average power and 
the peak power by a, 

where we assume < a < 1. Note that a = 1 corresponds to 
the case with only a peak-power constraint. Similarly, a<l 
corresponds to a dominant average-power constraint and only 
a very weak peak-power constraint. 

We denote the capacity of the described channel with 
peak-power constraint A and average-power constraint £ by 
C(A, £). The capacity is given by [3] 



C(A,£) = sup I(X;Y) 



(7) 



where I(X;Y) stands for the mutual information between the 
channel input X and the channel output Y, where conditional 
on the input x the output Y is Gaussian ~ A/"(a;,cr 2 ); and 
where the supremum is over all laws on X > satisfying 
Pr[X > A] = and E [X] <£. 

In the case of only an average-power constraint the capacity 
is denoted by C(£). It is given as in © except that the 
supremum is taken over all laws on X > satisfying E[X] < 
£. 

The derivation of the upper bounds is based on a technique 
introduced in [4]. There a dual expression of mutual informa- 
tion is used to show that for any channel law W(- 1 •) and for 
an arbitrary distribution R(-) over the channel output alphabet, 
the channel capacity is upper-bounded by 



C<E Q , [D{W{-\X)\\R{-))]. 



(8) 



Here, -D(-||-) stands for the relative entropy [5, Ch. 2], and 
Q*(-) denotes the capacity-achieving input distribution. For 
more details about this technique and for a proof of (O, see 
[4, Sec. V], [6, Ch. 2]. The challenge of using (|8]i lies in a 
clever choice of the arbitrary law R(-) that will lead to a good 
upper bound. Moreover, note that the bound ^ still contains 
an expectation over the (unknown) capacity-achieving input 
distribution Q*(-). To handle this expectation we will need to 
resort to some further bounding like, e.g., Jensen's inequality 
[5, Ch. 2.6]. 

The derivation of the firm lower bounds relies on the entropy 
power inequality [5, Th. 17.7.3]. 

The asymptotic results at high power follow directly by 
evaluation of the firm upper and lower bounds. For the low 
power regime we introduce an additional lower bound which 
does not rely on the entropy power inequality. This lower 
bound is obtained by choosing a binary input distribution, a 
choice which was inspired by [7], and by then evaluating the 
corresponding mutual information. For the cases involving a 
peak-power constraint we further resort to the results on the 
asymptotic expression of mutual information for weak signals 
in [8]. 

The results of this paper are partially based on the results 
in [9] and [6, Ch. 3]. 



The remainder is structured as follows. In the subsequent 
section we state our results and in Section [HI] we give a brief 
outline of some of the derivations. 

II. Results 

We start with an auxiliary lemma which is based on the 
symmetry of the channel law and the concavity of channel 
capacity in the input distribution. 

Lemma 1: Consider a peak-power constraint A and an 
average-power constraint £ such that a = ^ > i. Then the 
optimal input distributiorQ in (0 has an average power equal 
to half the peak power 

Eq*[A1 = 1 
A 2' 



(9) 



irrespective of a. I.e., the average-power constraint is inactive 
for all a E (f >l]> an d m particular 

C(A,aA) = c(A,^, ~<a<l. (10) 

We are now ready to state our results. We will distinguish 
between three different cases: in the first two cases we impose 
on the input both an average- and a peak-power constraint: in 
the first case the average-to-peak power ratio a is restricted 
to lie in (0, |), and in the second case it is restricted to lie 
in (Note that by Lemma Q] 4 < a < 1 represents 

the situation with an inactive average-power constraint.) In the 
third case we impose on the input an average power constraint 
only. 

In all three cases we present firm upper and lower bounds 
on the channel capacity. The difference of the upper and lower 
bounds tends to when the available average and peak power 
tend to infinity with their ratio held constant at a. Thus we can 
derive the asymptotic capacity at high power exactly. We also 
present the asymptotics of capacity at low power: for the cases 
involving a peak-power constraint we are able to state them 
exactly, and for the case of only an average-power constraint 
we give the asymptotics up to a constant factor. 

A. Bounds on Channel Capacity with both an Average- and a 
Peak-Power Constraint (0 < a < \) 

Theorem!: If < a < 4, then C(A,aA) is lower- 
bounded by 



C(A,aA)>-l 0g (^l + A ^ 2 {—- 
and upper-bounded by each of the two bounds 

C(A,aA)<ilog (l + a (l-a)^ 

C(A, aA) < (l - Q (^) - Q ( * +(1 ; a)A 

, [A e^-e-^+i) \ 
■ log — 



(ID 



(12) 



a V2^(l~ 2Q(f)) 



'it was shown in [1] that the optimal input distribution in Q is unique. 
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Fig. 1 . Bounds of Theorem [2] for a choice of the average-to-peak power 
ratio a = 0.1. The free parameters have been chosen as suggested in (16) 
and (T7J. The maximum gap between upper and lower bound is 0.72 nats (for 
A/tr w 11.8 dB). 
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Fig. 2. Bounds of Theorem [2] for a choice of the average-to-peak p ower 
ratio a = 0.4. The free parameters have been chosen as suggested in \16\ 
and (T7J. The maximum gap between upper and lower bound is 0.56 nats (for 
A/a « 7.1 dB). 
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Here Q (■) denotes the Q-function defined by 

/•°° 1 .2 

e(0 ^ __. e -Tdi, v^e 



/2tt 



(13) 



(14) 



j« > and 5 > are free parameters; and fi* is the unique 
solution of 

1 e-^ 
a = -. (15) 

H* 1 - e~» 

Note that /x* is well-defined as the function (i* i— > -=j- — 
7j-=7?t is strictly monotonically decreasing over (0, oo) and 



tends to ^ for it* | and to for //* | oo. 

A suboptimal but useful choice for the free parameters in 
the upper bound ( fl3l l is 

« = «(A)4 (7 logfl + ^ i (16) 

(17) 

where /i* is the solution to (fT3T >. For this choice and for a = 
0.1 and 0.4 the bounds of Theorem [2] are depicted in Figures Q] 
and |2 

Theorem 3: If a lies in (0, i), then 



X(«)= lim |c(A,aA)-log^| 



= -- log2^e - (1 - a)// - log(l - a M *) (18) 



and 



lim 



C(A,oA) _ a(l - a) 



(19) 



aIo A 2 /a 2 

B. Bounds on Channel Capacity with a Strong Peak-Power 
and Inactive Average-Power Constraint ( * < a < 1) 

By Lemma Q] we have that for i < a < 1 the average- 
power constraint is inactive and C(A, a A) — C(A, ^A). Thus 
we can obtain the results in this section by simply deriving 
bounds for the case a = |. 

Theorem 4: If a € [i , l], then C(A, aA) is lower-bounded 

by 

C(A,aA)> ibgfl+^-^V (20) 



2" ' 2irea 2 J ' 

and is upper-bounded by each of the two bounds 



C(A,aA)<ilog (l + A_ 



C(A,aA) < 1-2Q 



(21) 




A + 26 



cV^F (1-2Q(£)) 
-S?, (22) 



where S > is a free parameter. 
We suboptimally choose 

6 = 6(A) =a log 



(23) 



For this choice the bounds of Theorem |4] are depicted in 
Figure [3] 

Theorem 5: If a lies in [|, l], then 



and 



^ lim (c(A,aA)-log-l = 
C(A,aA) 1 



lim ■ 



-log27re (24) 



(25) 



A|o A 2 /ct 2 8' 
Note that (1241 and d25l > exhibit the well-known asymptotic 
behavior of the capacity of a Gaussian channel under a peak- 
power constraint only [3]. 
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Fig. 3. Bounds on the capacity of the free-space optical intensity channel with 
average- and peak-power constraints for a > ^ according to Theorem|4] This 
includes the case of only a peak-power constraint a = 1. The free parameter 
has been chosen as suggested in j23t . The maximum gap between upper and 
lower bound is 0.54 nats (for A/a ~ 7 dB). 
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Fig. 4. Bounds on the capacity of the free-space optical intensity channel 
with only an average-power constraint accor ding to T heorem [6] The free 
parameters have been chosen as sugg ested in (29) (32}. The maximum gap 
between upper and lower bound is 0.64 nats (for £ /a m 1.8 dB). 



C. Bounds on Channel Capacity with an Average-Power Con- 
straint 

Finally, we consider the case with an average-power con- 
straint only. 

Theorem 6: In the absence of a peak-power constraint the 
channel capacity C(£) is lower-bounded by 



C(«>ilo S (l + |^), 
and is upper-bounded by each of the bounds 



(26) 



C(£) < log pe~^ + V2iraQ 



log I v 2na 



2V27TCT 

1 



2a 2 



5 + e 



- - log 2irea , 



5>0, (28) 



where (3 > and S are free parameters. Note that bound ( l27l ) 
only holds for 5 < —oe~?-, while bound ( f28l only holds for 
8 > 0. 

A suboptimal but useful choice for the free parameters 
in bound ( |27| i is shown in d29l i and ( f30b and for the free 
parameters in bound (l28l i is shown in OTb and d32l at the top 
of the next page. For these choices, the bounds of Theorem|6] 
are depicted in Figure |4] 

Theorem 7: In the case of only an average-power con- 
straint, 

X£ = lim (c(£ ) - log -1 = \ log ± (33) 



and 



lim • 

£10 



C(£) 



c(g) 



<2, 



lim ■ 



> 



1 

71' 



(34) 



(35) 



Note that the asymptotic upper and lower bound at low SNR 
do not coincide in the sense that their ratio equals 2\/2 instead 
of 1. However, they exhibit similar behavior. 

III. Derivation 

In the following we will outline the derivations of the firm 
lower and upper bounds given in the previous section. 

One easily finds a lower bound on capacity by dropping 
the maximization and choosing an arbitrary input distribution 
Q(-) to compute the mutual information between input and 
output. However, in order to get a tight bound, this choice 
of Q(-) should yield a mutual information that is reasonably 
close to capacity. Such a choice is difficult to find and might 
make the evaluation of I(X;Y) intractable. The reason for 
this is that even for relatively "easy" distributions Q(-), the 
corresponding distribution on the channel output Y may be 
difficult to compute, let alone h(Y). We circumvent these 
problems by using the entropy power inequality [5, Th. 17.7.3] 
to lower-bound h(Y) by an expression that depends only 
on h(X). I.e., we "transfer" the problem of computing (or 
bounding) h(Y) to the input side of the channel, where it is 
much easier to choose an appropriate distribution that leads to 
a tight lower bound on channel capacity: 



C = sup I(X;Y) 
«(•) 
>I(X;Y) 



for a specific Q(-) 



(36) 
(37) 



8 = 8(S) = -2o-Wlog|, for - < e~A « -0.4 dB, (29) 

O (T 



l/_ (J \ 1 / / _ (J \ 2 / (J \ /— - J _ / 8 



" = m ^2{ £ + vrJ + 2n £ + vrJ +4 { £+ vrJ V2lT ^ Q U» • (30) 



* = ^<7log(l + -), (31) 



/= /Cf) 4 J. + £ + * e -j£^ + lJfs + g +° e -£*\ + ±fs + g + -?L e -£z ) v 2^r^G ( - ). (32) 



/i(A + Z)' 



> 



for a specific Q(- 
g 2h(X) , e 2h(Z) 



for a specific Q(-) 



for a specific Q(-) 



log 1 



e 2h(X) 

27recr 2 



(38) 
(39) 

W (40) 
(41) 



for a specific Q(-) 

where the inequality in (1401 follows from the entropy power 
inequality. To make this lower bound as tight as possible we 
will choose a distribution Q(-) that maximizes differential 
entropy under the given constraints [5, Ch. 12]. 

The derivation of the upper bounds in Section [TT] are based 
on the duality approach ((HJ. Hence, we need to specify a 
distribution R(-) and evaluate the relative entropy in (0. 

We have chosen output distributions R(-) with the following 
densities. For (fT2l we choose 



R'iv) 



i 



y/2n {^+£{k-£)) 
for ( fT3l we choose 



R'iv) = 



-e 2a 



V27TCT 

i_ M(i-2a(|)) 

A ' 



(»-A) 2 



/27T(T 



-g 2cr2+2£(A-£) ■ 

V < -6, 
-8 < y < A 
y > A + 8, 



(42) 



(43) 

where 8 > and /i > are free parameters; for (|2TT > we 
choose 



for ( l22l we choose 



1 



(-IT 



(44) 



2tt a 2 



27TCT 



-2Q(#) 

A+2<5 • 



y < 

-8 < y < A + 5, 



(45) 



where 8 > is a free parameter; and for d27l > and d28| i we 
choose 



R'(y) 



1 



"72- 



7 , y < -5, 



(46) 



/3e 2^+V27ro-Q(£) 

where <5 € K and /3 > are free parameters. In the derivation 
of (f28b we then restrict 8 to be nonnegative, while in the 
derivation of d27t we restrict <5 < — ere~2. 
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